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43. 


ON THE DIFFERENTIAL EQUATIONS WHICH OCCUR IN 
DYNAMICAL PROBLEMS. 


[From the Cambridge and Dublin Mathematical Journal, vol. 11. (1847), pp. 210—219.] 


JACOBI, in a very elaborate memoir, “Theoria novi multiplicatoris systemati æqua- 
tionum differentialium vulgarium applicandi”(‘), has demonstrated a remarkable property 
of an extensive class of differential equations, namely, that when all the integrals of 
the system except a single one are known, the remaining integral can always be deter- 
mined by a quadrature. Included in the class in question are, as Jacobi proceeds to 
show, the differential equations corresponding to any dynamical problem in which 
neither the forces nor the equations of condition involve the velocities; i.e. in all 
ordinary dynamical problems, when all the integrals but one are known, the remaining 
integral can be determined by quadratures. In the case where the forces and equations 
of condition are likewise independent of the time, it is immediately seen that the 
system may be transformed into a system in which the number of equations is less 
by unity than in the original one, and which does not involve the time, which may 
afterwards be determined by a quadrature*; and, Jacobi’s theorem applying to this new 
system, he arrives at the proposition “In any dynamical problem where the | forces 
and equations of condition contain only the coordinates of the different points of the 
system, when all the integrals but two are determined, the remaining integrals may 
be found by quadratures only.” In the following paper, which contains the demonstra- 
tions of these propositions, the analysis employed by Jacobi has been considerably 
varied in the details, but the leading features of it are preserved. 


1 Crelle, t. xxvm. [1844], pp. 199—268 and t. xxrx. [1845], pp. 213—279 and 333—376. Compare also the 
memoir in Liouville, t. x. [1845], pp. 887—346. 
2 For, representing the velocities by 2’, y’ ... the dynamical system takes the form 
E AE OAE A et 2 Oey wee AF es 


and the system in question is simply dz: dy...: dx’: dy’... =a :y'..: X: YY... 
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§ 1. Let the variables æ, y, z,... &c. be connected with the variables u, v, w, ... 
by the same number of equations, so that the variables of each set may be considered 
as functions of those of the other set. And assume 


da dy ...=Vdudv...; 


if from the functions which equated to zero express the relations between the two 
sets of variables we form two determinants, the former with the differential coefficients 
of these functions with respect to u, v,... and the latter with the differential coeffi- 
cients of the same functions with respect to æ, y,... the quotient with its sign changed 
obtained by dividing the first of these determinants by the second is, as is well known, 
the value of the function V. 


Putting for shortness 


da dy _ _ da _, dy _ 
a de ) , E A y ala dv E, &e 
du du dv , dv _ 

and da 7 A a ai da > i 


V is the reciprocal of the determinant formed with A, B,...; A’, B’,..., &c.; or it is the 
determinant formed with a, B,...@’, B’,..., &C. 


From the first of these forms, i.e. considering V as a function of A, B,... 
dV _ dV _ dV , TAN $ 
a dA’ =~ ¥%: dB’ 7 VP 


where the quantities a, B,...a’, B’,... and A, B,... A’, B’... may be interchanged pro- 
vided — V be substituted for V. (Demonstrations of these formule or of some equivalent 
to them will be found in Jacobi’s memoir “De determinantibus functionalibus,” Crelle, 
t. XXII. [(1841) pp. 319—359].) 


Hence 
dV +adA PPIE FEJA + RaR 4... =0, 


or reducing by 
dA dB , dA’ _dB' be 
Deia Geyer del’ a 
this becomes 


gave (S d+ dy + ..) +8 E dt dy +...) vee 


+a (E do + 2B ay + ie) +8 (e dy + Fe nee 
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or, reducing, 
dA dA’ dB dB oe i 
aN (au t Gy to) det (Gat Ge to) yt 203 
whence separating the differentials and replacing A, A’,...; B, B’,...; by their values 
1dV.ddu.d d 
Zde ude de de.de.. agp 
1dV dad. d ad 
va aay ay 
(in which —V, w, v...; æ, y... may be substituted for V, æ, y...; u, v...) 
§ 2. Let X, Y...be any functions of the variables a, y,... and assume 
du du 
U= Xx ih ae 
du dv , 
V= X Tt Eag reni 
U, V, ... being expressed in terms of u, v,.... Then 
dU LEN d du d dv d du d dv 
ee he =X (ie pto et) ¥ (apa taeag ts) 
dX du aX dv dY du dY dv 
' dU av dV T dX aY 
Also, whatever be the value of M, 4 
dMV dMV dMV adMV 
and from these two properties, 
adMV U aUVV 


dMX  adMY 
du * dy nev (Se A 


dx dy 
§ 3. Consider the system of differential equations 
de: dy: dz....=X:Y: 245 


6 


(where, for greater clearness, an additional letter z has been introduced). From these we 


deduce the equivalent system 


du:dv:dw...=U:V:W.... 
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Suppose that u and v continue to represent arbitrary functions of æ, y, z,... but that the 
remaining functions w, ... are such as to satisfy W=0,... (so that w,... may be 
considered as the constants introduced by obtaining all the integrals but one of the 
system of differential equations in 2, y, z, ...), we have 


dMVU | dMVV _ q /dMX , dMY , dMZ 
du ; mesa wt a TAN 


Also the only one of the transformed equations which remains to be integrated is 
du: dv=U: V, or Vdu— Udv=0, 


(in which it is supposed that U and V are expressed by means of the other integrals in 
terms of u and v). 


Suppose M can be so determined that 


dMX dMY dMZ 


Taal ape de T eaS, 


(M is what Jacobi terms the multiplier of the proposed system of differential equations) : 
then 
dMVU Ms dMVV _ 0 
du ae ie 
or MV is the multiplier of Vdu— Udv =0, so that 


fury (Vdu — Udv) = const. 


Hence the theorem :-—“ Given a multiplier of the system of equations 
apy: ay Gee ee AY 2... 


(the meaning of the term being defined as above), then if all the integrals but one of 
this system are known, the remaining integral depends upon a quadrature.” 


Jacobi proceeds to discuss a variety of different systems of equations in which it 
is possible to determine the multiplier M. Among the most important of these may 
be considered the system corresponding to the general problem of Dynamics, which 
may be discussed under three different forms. 


§ 4. Lagrange’s first form’. 


Let the whole series of coordinates, each of them multiplied by the square root 


of the corresponding mass, be represented by 2, y,... and in the same way the whole 
series of forces, each of them multiplied by the square root of the corresponding mass, 
by P, Q,...; then the equations of motion are 

do ` y- e 

a oe oA. 


1 I have slightly modified the form so as to avoid the introduction of the masses, and to allow «æ (for 
instance) to stand for any one of the coordinates of any of the points, instead of standing for a coordinate 
parallel to a particular axis. 
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where 
d® = 
d® dọ . 
Y=Q+ X Tu +y +..., 
T PRT dy 
where @=0, @=0,... are the equations of condition connecting the variables, and 


À, M+». coefficients to be determined by substituting the values of P’ &c. in the 


equations ot = 0, 


pendent of the velocities. 


=0, &c. It is supposed that as well P, Q,... as ©, ®,... are inde- 


In order to reduce these to an analogous form to that previously employed, we 
have only to write 


dx, dy 
eae Tak oe 


which gives 
dt: da: dy:dz...: da’: dy: dz’... 
et a Ad a SO ee ge ae 


Supposing that M is independent of a’, y, 2’,... the equation on: which it depends 
becomes immediately 


aM + M (Tet +a.. .) =0, 
d. 
where for shortness 
dt nk a Ee 
b=7t+e wt Gy +? da teers 
To reduce this we must first determine the values of à, w..., and for this we have 
LO dO dx dO dy 
de Teet D ye =0, &e. 
IF d® 
i.e. COFI a THH -+tarthwt+gv+...=0, 


® 
sp + pt Te + 95 +...thry+bu+fv +... =0, 


v 
ree Ete, +. -t+gry=fu+ov+...=0. 
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where for greater clearness an additional letter of the series ©, ®... has been introduced, 


and where 
s 2 
a= (T + e asa 


= (a) +(a,) Er 
A (2 BoB e 


Hence differentiating with respect to a’, 


sett ay i ae a. ao 


d® dà du dv 


dY 


dx du , 
yo a +Ig7ttagt Cay = 95 


ur representing by K the determinant formed with the quantities a, h, g,... h, b, fy... 
9, f, ¢,...and by A, H, G,... H, B, F,....G, F, C,... the inverse system of coefficients, 


we have 


d@ „dò „~ dY dn 
2 (43+ HBS + G3...) + KX =0, 
d® d „dY du 
2 (H + Bb T+ FBT...) + KM =0, 
d® „dð. „dY dv 
2 (G85 rete) =o, 
whence, multiplying by nA sA aa ... and adding, 
2 2 
as ($) +83 (T) naa rein, 


and, forming the similar equations with the remaining variables and adding, 


dX dY dZ )=0; 


dx’ tay? ae” 


dX dY aZ =o 


ie. èK + E(t art at 


Aba + Bob + Osc ... + 2Fèf + 2489 + 2Hèh +... + K( 
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Thus the equation in M reduces itself to 
KèM — MòK = 0, 


which is satisfied by M=K. It may be remarked that. K reduces itself to the sum 
of the squares of the different functional determinants formed with the differential 
coefficients of ®©, ®... with respect to the different combinations of as many variables 


out of the series a, y... 
§ 5. Lagrange’s second form. 


Here the equations of motion are assumed to be 


addaT aT 
di da’ — da —* ~° 
addaT dT 
dt dy dy ~ 
aa dT py 


where 2T represents the vis viva of the system, æ, y, Z,... are the independent variables 
on which the solution of the problem depends, and 2’, y, 2’,... their differential coeff- 
cients with respect to the time. It is assumed as. before P, Q, R... do not contain 


A i fs eae 
Suppose these equations give 
di: da:dy:dz...: de: dy : dz... 
mE ye aay ay Se ees 


then the equation which determines the multiplier M takes as before the form 


M+ (E +o tat .)=0, 
To reduce this equation, substituting for 7’ its value which is of the form 
T=} (av? + by? +02" ... + fyz + 2gza' + 2ha'y’ ...), 
and putting for shortness 
L=ax +hy +gz’..., 
M = ha! + by' + fe’ ..., 
N=ge'+fy +02’... 
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the equations which determine X, Y, Z... are 


aX +h¥ + 9Z...+8L -$ - P=0, 


AX + bY + ff... +M- To- R=, 


9X +fY +0L...+3N-T—R=0. 
Hence, differentiating with respect to 2’, 
aX ,dY, dZ > 
a gg teag t I gg + 8a =e (), 
aX  ,dY. 1.dZ dM dL _ 


dX dY dZ dN dL 
Bae hs ep hia: ar am 0; 


or representing by K the determinant formed with a, h, g,...h, b, f,...9, f, c,... and 
by A, H, G,... H, B, F,... G, F, © ... the inverse system of coefficients, we have 


dX F dM oah (dN an _ 
K & + Ada + Hoh + Gog... + af) a(g g). 
and similarly 
dY dL dM 7 dN dM _ 
K gy + Hoh + Bb + Fe... +H (Fy - o) + +F(G -3)-=0 
dZ dL dN\, „(dM dN 
K g + Gog + Ff + be... + G (T Tr) EGAR) i ci 


Hence, adding, 


dX dY dZ ; TEP 
K (qa + qy + qz =) + 43a + BOb + O8e ... +2F8/ + 208g + 2HBh ... = 0: 


and thus we have as before, though with symbols bearing an entirely different signification, 


dX dY aZ 
EGA T n.) + 8K =0; 


and thence K8M— MSK =0, and M = K. 
36—2 
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(The value of K in this section may I think be conveniently termed “the deter- 
minant of the vis viva,’ with respect to the variables a, y, z,.... It may be remarked 
that “the determinant of the vis viva” with respect to any other system of variables 
u, v, w, ... is =V?K, V as before.) 


§ 6. Third form of the equations of motion. [Hamiltonian form.] 


Here writing 
dT 
dat -E3 iy nee 


and taking t, a, y,...& 7,... for the variables of the problem [and considering T to be 
expressed as a function of these variables: to denote this change it would have been 
proper to use instead of T a new letter H] the equations of motion reduce them- 


selves to 
a. ae ON ee 


dn _ at dy dT. 
Ber get. eas 


E _aT p da_ar 


or putting for shortness 


dx dé 

aT dT 
r Ta dn ? 
i ; 


they become 
dt : da: dy: dz.,.: dE: dy: det 


re we et eee A 


Hence writing the equation in M under the form 


dE dH ax dY 
èM OF cc ee ar Jai 
+M (T+ eter at ) 0: 

d d He. a re 


we see immediately that (P, Q ... being as before independent of the velocities, and con- 
sequently of & n, €,...), 


ss ok -0, T+ dY 


+ aE Sane Se 


Hence 8M = 0, which is satisfied by M = 
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